Random changes of flow topology in two dimensional and geophysical turbulence 
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We study the two dimensional (2D) stochastic Navier Stokes (SNS) equations in the inertial limit 
of weak forcing and dissipation. The stationary measure is concentrated close to steady solutions of 
the 2D Euler equation. For such inertial flows, we prove that bifurcations in the flow topology occur 
either by changing the domain shape, the nonlinearity of the vorticity-stream function relation, 
or the energy. Associated to this, we observe in SNS bistable behavior with random changes from 
dipoles to unidirectional flows. The theoretical explanation being very general, we infer the existence 
of similar phenomena in experiments and in models of geophysical flows. 
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The largest scales of turbulent flows are at the heart of 
a number of geophysical processes : climate, meteorology, 
ocean dynamics, the Earth magnetic field. The Earth 
is affected on a very large range of time scales, up to 
millennia, by the structure and variability of these flows. 
Many of these undergo extreme and abrupt qualitative 
changes, seemingly randomly, after very long period of 
apparent stability. This occurs for instance for magnetic 
field reversal for the Earth or in MHD experiments [H, 
for 3D flows [2|, for multiple equilibria of atmospheric 
flows [K], for 2D turbulence experiments 0, 01 and for 
the paths of the Kuroshio and Gulf Stream currents @. 

Understanding these phenomena requires a statistical 
description of the largest scales of turbulent flows. Very 
few theoretical approaches exists due to the prohibitively 
huge number of degrees of freedom involved. Fruitful 
hints may be drawn from qualitative analogies with bista- 
bility in system with few degrees of freedom perturbed 
by noise [7|. However the range of validity of this ap- 
proach remains a tough scientific issue, because of the 
complexity of turbulent flows. What is the good theoret- 
ical framework for such phenomena ? In the following, 
we argue that 2D turbulence, because of its relative theo- 
retical simplicity, is a very interesting framework in order 
to address such an issue. 

In this letter we predict and prove the existence of ran- 
dom switches from dipoles to unidirectional flows (see Fig 
[1]), in the 2D Navier Stokes Eq. with random force (SNS). 
Similar random changes have already been observed in 
rotating tank experiments for quasi-geostrophic dynam- 
ics [31] . Following analogous theoretical considerations as 
for SNS Eq., we infer that such changes will generically 
occur within a large class of models like quasi-geostrophic 
(QG) or shallow-water (SW) models that describe atmo- 
spheric [K], and oceanic 0| large scales. The recipe we 
propose is to exhibit bifurcation lines representing abrupt 
change in steady solutions in the inertial limit and then 
look for the corresponding transitions in real flows. 

Geophysical and 2D inviscid flows are characterized 
by the conservation of energy and an infinite number of 
quantities (Casimirs), such as enstrophy. This property 
prevents direct energy cascade towards the small scales, 
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Figure 1: Time series and probability density functions 
(PDFs) for the order parameter z\ (see page 4) illustrating 
random changes between dipoles and unidirectional flows. 



by contrast with 3D turbulence. Then, the first phe- 
nomenon is an inverse energy cascade towards the large 
scales and a direct enstrophy cascade. Kraichnan classi- 
cal theory (gj studies the self-similar processes associated 
with these two cascades (see the recent spectacular dis- 
covery of conformal invariance consequences for the in- 
verse cascade [9]). The second phenomenon, the self orga- 
nization of the flow into jets and vortices, occurs if energy 
is not dissipated before reaching the largest scale. Then 
coherent structures break the self-similarity so that their 
study cannot be properly addressed using Kraichnan the- 
ory [Sj. A second classical theory, the so-called Robert- 
Sommeria-Miller (RSM) equilibrium statistical mechan- 
ics [lOj , predicts the self-organized structures for invis- 
cid decaying turbulence. However, this inviscid theory 
does not take into account the long-term effects of forc- 
ing and dissipation as well as the slow dynamics of the 
flow. Therefore, random changes of flow topologies can- 
not be explained by these two classical theories. 

As an alternative theoretical approach, we study sta- 
tistically stationary states of SNS Eq. Note thata self- 
similar growth of a dipole has been studied in em- 
phasizing transient growths: both approaches comple- 
ment each other. SNS Eq. on a doubly-periodic domain 
V = (0; 2tt5) x (0; 2tt) reads 



duo 



+v«Vw = — aoj+isAoj+yfar], v = e 2 xV^, uj = Aip, (1) 



where u, v and ip are respectively the vorticity, ve- 
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locity and stream function ; a is the Rayleigh fric- 
tion coefficient and v the viscosity. The force curl is 
r] = 7 /k^k (i) exp(zk.x)/ (27r), with {rj^} independent 

— 'k 

Gaussian white noises : (£')) = ^kk /( ^(^ — 

We impose £?o = ^ k |/k| 2 / |k| 2 = 1 so that a is the 
average energy injection rate. 

Euler Eq. (a = v = a = 0) conserve the kinetic energy 
E and vorticity moments Q n (Q2 is the enstrophy) 

E = - [ d 2 xv 2 and fi n = / d 2 xuo n . (2) 

Application of Ito formula to the energy, and averaging 
over the noise, leads to d (E) /dt= —2a (E) + a — v (^2). 
If (.) s denotes averages over the stationary measure, we 
have 2a (E) s + v (^2)5 = o~- It expresses the balance be- 
tween energy injection and energy dissipation. Clearly, 
for flows with energetic large scales, Rayleigh friction 
dominates dissipation 2a (E) s ^> v (^2) $ (mathemati- 
cally we consider the limit v — > for fixed a and assume 
^(^2)5 — > 0). It is natural to fix the average energy 
to be of order 1 by using a typical turnover time as a 
new time unit. Putting t' = y 1(7 1 (2a)t, uo 1 = \j2aj 'cruo, 

a'=(2af /2 /(2a 1 / 2 ) and v' = v(2a/a) 1/2 and dropping 
the primes, the dimensionless Eq. are 

^ + v «V uo — —auo + is Auo + V2ar]. (3) 
at 

The energy balance now reads (E) s + (is /2a) (^2)5 = 1- 
In these dimensionless unit the Reynolds number is 1/is 
and the Rayleigh number is R a = O (v -V uo / 'auo) = 1/a 
(arresting the inverse cascade before energy reaches the 
largest scale would requires a > 1). For most geophysical 
flows and experiments the case of weak forcing and dissi- 
pation is the most relevant one. We thus study the iner- 
tial limit «<1 (more precisely the limit lim^^o lim^-^o)- 
Without Rayleigh friction (a = 0), the previous dis- 
cussion is meaningless and the balance relation becomes 
2is (^2)5 = cr. By a natural time unit change, we can fix 
(Q 2 )s = 1. The nondimensional equation is then 

^+v-Vw = v Auj + V2urj. (4) 
at 

From a physical point of view, this last model is less 
relevant than ([3]) but is still very interesting from an 
academic point of view. A series of recent works has 
proved the existence of invariant measures, validity of 
the law of large numbers, central limit theorems, ergod- 
icity and some properties of stationary measures in the 
inertial limit v — > 0, balance relations (see [12] and refer- 
ences therein). All following considerations are relevant 
for both models (|3|4|) . in their respective inertial limits. 

We know since decades from real [13] or numerical 
[!, Q3] experiments, that for times large compared to 
the turnover time but small compared to the dissipation 
time, the largest scales of 2D Navier-Stokes turbulent 
flows converge towards steady solutions of Euler Eq. : 

v • V uj = or equivalently uo = / (ip) . (5) 



It appears to be true as well for the Euler Eq. For in- 
stance, RSM theory predicts / from given initial condi- 
tions. Given this empirical evidence, it is thus extremely 
natural to expect that in the inertial limit, measures for 
SNS are concentrated near steady Euler flows. We show 
numerical evidences of this fact in the following. 

The ensemble of steady Euler flows is huge, as it is 
parametrized by the function /. It will be proven that 
when either / or the domain shape is changed, bifurca- 
tions may occur. Such abrupt transitions lead to strong 
qualitative changes in the flow topology. In this criti- 
cal regime and under the action of a small random force 
in SNS, the system switches randomly from one type of 
topology to another. In the following, we show that this 
scenario is valid. 

We study a bifurcation diagram for stable steady Euler 
solutions, by considering 

S(E) = sup{S[uo] = [ d 2 x s(uo) \ S(uo) = E}, (6) 

uj JV 

where S{E) is the equilibrium entropy, S the entropy of 
uo ; the specific entropy s (uo) = — uo 2 /2-\~Yl n>2^2n/2 n uo 2n 
is concave assuming s even for simplicity. Critical points 
of (|6]) verify uo = f(i/j) = (s 7 ) -1 (— where f3 is the 
Lagrange multiplier associated with energy conservation. 
They are thus steady Euler flows, satisfying ([5]), and the 
knowledge of / or s are equivalent. From Arnold's the- 
orems [15| or its generalization, maxima of ([6]) are dy- 
namically stable. One can also prove that any solutions 
for ([6]) are RSM equilibria [16]. Even if it seems appeal- 
ing, there are no clear theoretical arguments for giving 
a thermodynamical interpretation to (|6]) in the SNS out- 
of-equilibrium context. We thus consider ([6]) only as a 
practical way to describe ensembles of stable steady Eu- 
ler solutions. 

Dipoles and unidirectional flows ("bars") have been ob- 
tained numerically Jl7|| as entropy maxima for 2D Euler 
Eq. with periodic boundary conditions, assuming sinh, 
tanh and 3-level Poisson uo — ip §5§ relations. According 
to [17| "which has the greater entropy (between dipoles 
and bars) depends on seemingly arbitrary choices". 

The fact that both unidirectional flows and dipole may 
be equilibria can be understood from the small energy 
limit of (0. Let us call {ei}i>i the orthonormal family of 
eigenfunctions of the Laplacian — Ae^ = A^, {eiej) v — 
Sij (< . >t>= j v dx and are arranged in increasing 
order). We decompose the vorticity as uo = X^i 6 ^^- 
The energy is then 2£{uo) — J2i>i \ luj 1- Since £{uo) is 
always positive (oo 2 ) v is small in the limit E — > 0, and 
only the quadratic part of S [uo] is relevant. Long but 
straightforward computation of (0 in the limit E — > 
gives 

uj - (2X 1 E) 1/2 e 1 with S(E) = -X 1 E + 0(a 4 XlE 2 ). (7) 

E — »0 

We thus conclude that the eigenmode with the small- 
est eigenvalue is selected, corresponding to the heuris- 
tic idea that energy condensate to the largest scale. 



For instance when the aspect ratio 5 > 1, the mode 
e\ = rii sin [(x + fa)/8\ is selected. This corresponds to 
a unidirectional velocity field vi = ri\ cos \{x + fa)/S] e y 
where fa is a phase associated to the translational in- 
variance. A dipole is actually a mixed state ae\ + (3e2 
with e2 = ri2 sin(y + fa) • In the weak energy limit, it can 
be selected only for the degenerate case Ai = A2. This 
happens for the square box 5 = 1. In such a case we 
can prove that the degeneracy is removed by the contri- 
bution of higher order terms in ([7j). From ([7j), we con- 
clude that the domain shape (ei) selects the equilibria 
for A2 — Ai ^> a^XfE whereas for A2 — Ai <C a^XfE the 
degeneracy is removed by the nonlinearity of / (a^). In 
order to study the bifurcation between these two behav- 
iors, we define g by A2 — Ai = gE and we study the 
small energy limit E — > 0, with fixed g. Straightforward 
computations lead to 

S(E) = -X ± E + E 2 max h(X), (8) 

with h{X) = {e\) v a 4 \j -gX + 2 7 A?a 4 X(l - X), where 
7 = ^{e\e2) v — ( e i) v > 0- The vorticity equilibria 
is then uo eq ~e^o \/2\ 1 E{l - X M )ei + \J2E\\X M ^2 
where Xm is the maximizer of h in (|S|). For Im = 
or Xm = 1, oj eq is an unidirectional flow whereas for 
< Xm < 1 it is a dipole (symmetric for Xm = 1/2). 
The selection occurs via maximization of h. When max- 
imizing h, the sign of the parameter a 4 plays a crucial 
role. We note that a 4 is intimately related to the shape 
of the relationship uj = f(ip) = (V) -1 {—f3ijj) t Indeed 
(V) -1 (—x) = x + a^x 3 + o(x 3 ) and when 0,4 > (resp. 
a 4 < 0), the curve /(^) bends upward (resp. downward) 
for positive i/j similarly to sinh (resp. tanh). 

The bifurcation diagram is summarized in Fig. [2] a). 
In the degenerate case (g = 0), the dipole is selected 
for <24 > (sinh like), whereas unidirectional flows are 
selected for a 4 < (tanh like). The term — gX favors 
the pure state e\ (X = 0). For a 4 < the unidirectional 
flow ei is always selected. More interestingly, for C14 > 
a bifurcation occurs along the critical line g* = 27A^<24 
between dipole and unidirectional flows. 

We have obtained the bifurcation diagram in the limit 
of small energy using the scaling A2 — Ai = gE. From a 
practical point of view, it is more convenient to work for 
a fixed aspect ratio 6. Using the relation g = ( A2 — Ai ) /E, 
we obtain that the critical line in a E — a 4 plane is the hy- 
perbola a^E = 8tt 2 (6 — 1) /3+o(5— 1). We use a continu- 
ation algorithm in order to numerically compute solution 
to (0 corresponding to f aA { x ) ~ (1/3 — 2a4) tanhx + 
(2/3 + 2*24) sinh x. The inset of Fig [2] b) shows good 
agreement for transition lines obtained either with the 
continuation algorithm or the low-energy limit theoreti- 
cal result, for S = 1.01. Figure[2]b) shows the bifurcation 
diagram for 6 = 1.1; in such a case the transition line is 
still very close to an hyperbola provided energy is small. 




Figure 2: Bifurcation diagrams for steady Euler flows a) in 
the g-d4 plane b) obtained numerically in an E — cu plane 
with S = 1.1. The inset illustrates good agreement between 
numerical and theoretical results in the low energy limit. 

Following the same reasoning, small-energy bifurca- 
tion diagrams could be computed for any Euler-like 
model like QG or SW models. Most often, the do- 
main shape selects the flow topology. When domain 
shape is varied, we meet eigenvalues degeneracy. In 
all these cases a bifurcation diagram can be computed 
where the transition line corresponds to the competition 
between the uj — ip ([5} nonlinearity and the domain shape. 

We expect to observe both dipoles and unidirectional 
flows in SNS. Numerical simulations in a square domain 
6=1 exhibit statistically stationary uo with a dipole 
structure (Fig. [3] a)), whereas for 6 > 1.1, nearly uni- 
directional flows are observed (Fig. [3]b)). This result 
has been confirmed both for a = and «/0, and for 
different v values and forcing spectra. One observes in 
Fig. [3] a uj — ip relation qualitatively similar to a sinh, in 
the dipole case and to a tanh in the unidirectional case. 
This confirms that uj remains close to steady Euler flows. 
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Figure 3: uj — ip scatter-plots (cyan). In black the same after 
time averaging (averaging windows 1 < r < 1/^, the drift 
due to translational invariance has been removed) a) dipole 
case with S = 1.03 b) unidirectional case 5 = 1.10. 

A very natural order parameter is \zi\, where z\ = 
(uj(x : y) ex.p(iy)) v . Indeed, for unidirectional flow 

uj = aei, z\ = 0, whereas for a dipole uj = a {e\ + e<i)^ 
\z\\ = a. Fig. [1] shows \z±\ time series for S = 1.02 and 
6 = 1.04. The remarkable observation is the bimodal be- 
havior in this transition range. The switches from \z±\ 
values close to zero to values of order of 0.6 correspond 
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to genuine transitions between unidirectional and dipole 
flows. The PDF of the complex variable z\ (Fig. [T]) ex- 
hibits a circle corresponding to the dipole state (a slow 
dipole random translation corresponds into to a phase 
drift for zi, explaining the circular symetry). The zonal 
state corresponds to the central peak. As S increases, 
one observes less occurrences of the dipole. For larger 
(resp. smaller) values of S only unidirectional (dipole) 
flows exist. The transition is also visible in other physical 
variables. For instance ^4 = ( la;] 4 )^ switches between a 
state with weak variance and low mean value (unidirec- 
tional) to an intermittent state with large variance and 
larger mean value (dipole). Topology changes are very 
slow dynamical processes : for the model ((4]), an aver- 
age transition time is of order \jv. For instance Fig. Q] 
represents 3.10 4 turnover times. For this reason, because 
of numerical limitations it has not been yet possible to 
obtain convincing analysis of the switch time statistics. 

In the spirit of [7] we look for low-dimensional analo- 
gies. When looking how evolves the PDF for the order 
parameter, while changing the control parameter, the 
dipole-unidirectional transition has striking similarities 
with the stochastic differential equation 



dx = x(n + x 2 — x 4 )dt + adW. 



(9) 



The deterministic part of (0 is the normal form for a 
generalized subcritical pitchfork bifurcation. For \i < 
— 1/4, one has a single stable fixed point x* = 0. For 
\i > there are three fixed point, one unstable xo = and 
two stable x h2 = ±(l + (l+4/i) 1 / 2 ) 1 / 2 . For fi e] -1/4, 0[, 
three stable fixed points coexist (xo and #1,2) and two 
unstable ones x 3A = ±(1 - (l + 4/i) 1/2 ) 1/2 . With additive 
noise (a 7^ 0), when /i < —1/4, the \x\ PDF has a single 
peak centered at x = 0. In the interval corresponding 
to /i e] — 1/4, 0[, an additional peak appears related to 
|^i,2 1- Finally, there is a transition for fi larger than 
and only one peak corresponding to \x\\ remains. 

However, we stress that a low-dimensional model like 
(0, as useful as it may be, lacks part of the phenomena. 
For instance, it can not explain why ^4 is intermittent 
while \z\ \ is not. Moreover, the role of turbulence here is 
not only to act as noise, but also to build up the large- 
scale flow by inverse cascade. The inverse cascade proper- 
ties are strongly affected by the existing large-scale flow, 
leading to the observed self-organization process. From 
a theoretical point of view, the main issue, beyond the 
scope of this letter, is to explain which of the Euler steady 
states will be selected by turbulence and to predict the 
relative frequency of such states. We thus need an alter- 
native theoretical approach bridging the gap between the 
two classical theories : self-similar inverse energy cascade 
on one hand and RSM equilibrium statistical mechanics 
on the other hand. 



In this letter, we have not addressed the other cru- 
cial issue : fluctuations. Some very interesting results 
and considerations on small-scale fluctuations for turbu- 
lence dominated by large-scale flows may be found in 
jllL [Til . fl9| . In forthcoming works, the statistical prop- 
erties of these random change of flow topologies and of 
fluctuations will be investigated. Finally, it will also be 
extremely interesting to analyze the connexions with sim- 
ilar transitions observed in other contexts 

We infer similar random flow topology changes for 
other geometry for 2D SNS, QG and SW models. Us- 
ing simple generalization of our analysis, rotating tanks 
experiments can be designed in order to observe similar 
phenomena. This study also suggests that flows like the 
Kuroshio currents [6] or the Gulf Stream might be close 
to steady solutions of inertial models. 
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